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We investigate singly and doubly charged flux tubes in U(l) lattice gauge theory By 
' simulating the dually transformed path integral we are able to consider large flux tube 

f — , lengths, low temperatures, and multiply charged systems without loss of numerical precision. 

We simulate flux tubes between static sources as well as periodically closed flux tubes, 
calculating flux tube profiles, the total field energy and the free energy. Our main results are 
that the string tension in both three and four dimensions scales proportionally to the charge 
^ - which is in contrast to previous lattice results - and that in four-dimensional U(l) there 

Q\ ■ is an attractive interaction between flux tubes for (3 approaching the phase transition. 



1 Introduction 

The formation of colour-electric flux tubes provides an intuitive physical picture for perma- 
nent confinement of quarks within hadrons. Moreover, this picture seems to be realized within 
QCD, as has widely been observed in lattice simulations. The mechanism leading to flux tube 
formation, however, is still subject to analytical and numerical investigations. A wide class 
of phenomenological models incorporates the existence of flux tubes: bag models pj, models 
describing the QCD vacuum as a perfect dia-electric medium |^], and the physically appealing 
picture of a dual superconductor ||]. The question arises which of these effective models are 
compatible with the numerical results of lattice simulations. 

The shortcoming of such lattice calculations, however, is the fact that regarding sufficiently 
long flux tubes, i. e. large lattice sizes, becomes a numerically very difficult task. Thus many in- 
teresting questions could not be answered until now with sufficient precision. Numerical troubles 
become even more serious, if one considers sources in higher representations which is essential 
to discriminate among the various phenomenological models 0]. Because of these difficulties it 
is convenient to consider gauge groups simpler than SU(3), like SU(2), or even compact QED 
(U(l) gauge theory) which also exhibits confinement of charges and can be used as a prototype 
for investigating the formation of flux tubes. 

In ref. Q SU(2) and U(l) lattice gauge theory in three dimensions were used to calculate 
some features of flux tubes generated by sources in different representations. In particular, the 
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cross section of flux tubes, the string tension and the spatial distribution of the energy density 
were investigated. The result was obtained that the cross section A is nearly independent of the 
representation j, and the string tension is proportional to the eigenvalue of the Casimir operator, 
i. e. a oc j(j + 1), which corresponds to the squared colour charge Q| (this has already been 
observed in previous lattice calculations). For multiple charges Q = ne in three-dimensional 
U(l) it was found that a oc n 2 . The authors conjectured that their results should also be relevant 
for flux tubes in four-dimensional QCD; the similarity of results obtained in SU{2) and U(l) 
was interpreted as indirect support for the dual superconductor picture of confinement, which 
is well established for U(l) lattice gauge theory. 

However, as was also realized in ref. [gj], this situation is in conflict with some of the above 
mentioned phenomenological models. According to our opinion, even the dual superconductor 
scenario is not compatible with the U(l) result a oc n 2 of ref. Q. We expect that in a dual 
superconductor the string tension increases linearly with the number of flux quanta, apart from 
interaction between flux tubes: In a type-II superconductor a doubly charged flux tube will split 
into two separated flux tubes. The energy will be just twice the energy of a single flux tube, if 
their transverse distance from each other is sufficiently large. In a type-I superconductor there 
is a negative interaction energy for two flux tubes which should diminish the ratio of string 
tensions below two. 

In fact we will demonstrate that the string tension scales proportionally to the charge rather 
than the squared charge in the confinement phase of U(l). In addition, investigating the inter- 
action energy of doubly charged systems we will discuss the question if four-dimensional U{\) 
behaves like a type-I or a type-II superconductor. The results concerning the flux tube cross- 
section and the energy obtained in ref. [|j on the other hand remind of the classical Coulomb 
case: The field strength is doubled everywhere in space, if we replace single by double charges, 
and the energy scales like the squared charge. 

It is the main aim of this paper to test the behaviour of flux tubes with respect to the 
above mentioned attributes within the dual formulation of four-dimensional U (1) lattice gauge 
theory. Simulating the dually transformed theory |J provides two great advantages compared 
to standard lattice simulations: The obtained signal does not decrease with the charge distance 
(the length of the flux tube), and a simulation with multiply charged sources can be performed 
with equal accuracy. Thanks to these features we will be able to perform a detailed investigation 
of single and double flux tubes in four-dimensional U(l), and also analyze the reasons for the 
discrepancy of our results to those in ref. 0] for the three-dimensional case. 

The remaining part of this article is organized as follows: In section 2 we shortly review the 
techniques used for the simulations. They are based on the method introduced in ref. Q, in 
addition we explain how to implement periodically closed flux tubes and doubly charged sources, 
and how to calculate the free energy. Then, in section 3, we present our results for flux tubes 
between static (double) charges in four-dimensional U(l). In section 4 we calculate the energy 
of periodically closed flux tubes (torelons). After a detailed analysis of the three-dimensional 
case in section 5, we finally turn to an investigation of the interaction between flux tubes in the 
four-dimensional theory in section 6. 
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2 Methods 



We are interested in the expectation values of physical observables O in the presence of a static 
charge pair at distance d, which are determined by 



( L(0)L+(d)O(x )) 
(L(0)L+(d)) 



where L denotes the Polyakov loop L(r) = rjfc=i ^x=(f,fca),/i=4 an d the angle brackets mean the 
averaging according to the path integral using the standard Wilson action 

S W = I3 [l-coed0]. (2) 

plaquettes 

(5 = 1/e 2 is the inverse coupling, and d6 is the discretized exterior derivative of the phases 8 
of the link variables U x „ = e lSx ^ and is assigned to a plaquette of size a 2 . For convenience we 
choose the notation of lattice differential forms. 

Each of the three expectation values appearing on the right hand side of (|l]) can be dually 
transformed as discussed in ref. M. The path integral including two external charges (£+ and 
£_ denote the world-lines of a static charge pair) reads 

Z QQ = U r M exp^^-cos^n^IF 9 = W iN TU2Tl e ~ Pl \\d*i + *n\m, 

C+ C- *l 

(3) 

where the first product on the right hand side is performed over all links of the dual lattice, and 
the Boltzmann factor is a product of modified Bessel functions over the dual plaquettes. *l is 
an integer valued 1-form. It is kept as dynamical variable and updated in simulations, while the 
2-form *n is fixed to be equal to 1 on the minimal straight area between the two Polyakov loops, 
and zero everywhere else. If one interprets the dually transformed partition function as a certain 
limit of a dual Higgs model ||, n has the meaning of a dual Dirac sheet, and *k = d*l + *n is 
the field strength. For the U(l) path integral without external sources one obtains the same as 
expr. (||), but with *n = 0. 

Further, we perform the duality transformation for the numerator in the correlation function 
(|l]). For the determination of the flux Fn through a given plaquette Do, we have to insert the 
operator Fn Q = y/]3 sin d9n || 0] into the path integral in addition to the Polyakov loops. One 
obtains 

/ \ / 1 \ dual 

( F °o) Q Q = Z QQ ^Ue- i9 Ue ie ^ S \nd9a J = \^* k °°) QQ ' ^ 

where the angle brackets in the final expression denote averaging over the dual degrees of freedom 
according to partition function (|3]), and *k = d*l + *n. The vacuum expectation value of the field 
strength F is of course zero. For the square of the field strength, which we need for determining 
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the local energy and action density, we take the usual operator /2 = (5 cos d#n []. In the dual 
formulation this becomes 

I (F 2 v _ / \" _ / MS. \" ( 5) 

An identity equivalent to (|5|) was already derived in ref. || for calculating the electromagnetic 
energy density for flux tubes in three-dimensional U(l) lattice gauge theory, using the polymer 
formulation of the path integral. 

As pointed out in ref. @, the correspondence of the dually transformed U(l) path integral 
to a limit of a dual Higgs model can qualitatively explain the formation of a flux tube (as 
product of a fluctuating string), but it does not help very much for interpreting the U(l) data in 
terms of a classical dual superconductor. Nevertheless, the above introduced dual formulation 
of expectation values in the presence of external charges provides a very efficient tool for flux 
tube simulations for several reasons. 

Contrary to standard lattice simulations, the confinement phase is the weakly coupled one in 
the dual theory, therefore there are less quantum fluctuations. Even more important is the fact 
that it is not necessary to project the charge-anticharge state out of the vacuum: On a lattice of 
same size charge pairs with arbitrary distance can be simulated with equal accuracy. For each 
charge distance d a separate simulation has to be performed, and also for the subtraction of the 
vacuum expectation value in (|l]) resp. (||) a distinct run is needed. The main advantage for the 
aims of this paper lies in the implementation of doubly charged flux tubes by just setting the 
integer field n on the dual Dirac sheet to 2 instead of 1. 

Our typical measurements are taken from 5 • 10 5 — 2 • 10 6 configurations (depending on 
the lattice size), where each 10th configuration is evaluated. The updating procedure uses a 
standard Metropolis algorithm; for determining the errorbars we took blocks of 100 evaluated 
configurations. We investigated the whole range of inverse coupling (3 in the four-dimensional 
theory. To compare the results in the confinement phase to those in the Coulomb phase, we 
focus on simulations at f3 = 0.96 and (3 = 1.10. In three-dimensional lattice U(l), where there 
is no phase transition to a deconfined phase, we also analyzed the /3-dependence of the string 
tension. 

For the required spatial size of the lattice there is the simple requirement that its spatial 
extent has to be greater than the extent of the considered flux tube both in longitudinal and in 
transverse direction in order to avoid finite size effects. For our simulations we used a lattice size 
of 8 3 x 32 for charge distances d = la — 5a, 12 3 x 32 for simulating d = 5a — 8a, and 16 3 x 32 for 
d = 8a — 12a. Further we performed some simulations on larger lattices; the largest investigated 
charge distance was d = 22a (on a 26 3 x 32 lattice). Finite temperature effects may be neglected 
if the temporal extent of the lattice exceeds the length of the flux tube (note that enlarging the 
temporal extent does not lead to a fall-off in the signal) ||]. A comparison between results from 
standard and from dual simulations has been performed in ref. ||. 

1 In Minkowski space the electric energy density is then given by = the magnetic energy density 

by Bf/2 = -FI.J2. 
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It turns out to be very convenient for our purposes to study flux tubes not only between 
static charges but also wrapping around the periodic lattice. These closed flux tubes ("torelons") 
are free of end effects and allow for a more reliable calculation of the string tension as well as 
the contribution from string fluctuations [|n]]. In standard simulations such torelons can be 
implemented by spatial Polyakov loops which are sufficiently separated in time direction. In 
the dual formulation it is straightforward to "remove the sources": The sheet defined by n = 1 
(resp. 2 for a double flux tube) is winding around the lattice both in temporal and longitudinal 
direction, and 5n = everywhere. In such a simulation the flux tube is not fixed any more and 
may move in transverse directions. Therefore we do not calculate field profiles, but just measure 
the energy of this state compared to the vacuum. 

Finally we have to specify what we mean by "energy". One approach is to sum up the 
spatial energy density, calculated from the expectation values of the squared fields. This is a 
very difficult task in standard lattice simulations, but within the dual simulation we are able 
to determine the total field energy as a function of the flux tube length. On the other hand, 
the usual "potential" is the zero temperature limit of the free energy relative to the vacuum, 
calculated from the correlation function between two Polyakov loops. Because this corresponds 
to a ratio of dual partition functions, and not to a dual expectation value, it cannot be obtained 
directly by a dual simulation. Nevertheless, we may relate the derivative of the free energy with 
respect to (3 to the average total action of the system: 

(6) 

Therefore it is possible to determine the free energy by a numerical integration over the average 
action relative to the vacuum as a function of /3, starting from the well-known strong coupling 
behaviour Fq^q^) ~* ~d ln(/5/2) as (3 — > 0. In the three-dimensional theory this procedure 
yields reliable results for the whole /3-range, whereas in four-dimensional U(l) integrating over 
the peak in the action density at the phase transition point (/? ~ 1.01) leads to large errorbars 
in the Coulomb phase. Above the phase transition it is therefore preferable to perform the 
integration starting from the weak coupling limit of U(l) gauge theory. It is obvious that the 
free energy for periodically closed flux tubes can be calculated in the same way. 

We may recapitulate that by dual simulations we are able to investigate large charge distances 
(flux tube lengths) at very low temperatures. Simulating doubly charged systems is possible 
without any loss of accuracy, and by studying periodically closed flux tubes one avoids being 
disturbed by end effects. For our analysis we will consider both the total field energy and the 
free energy of flux tubes in U(l) lattice gauge theory. 

3 Flux tubes between static double charges 

Usually a Wilson loop or a pair of Polyakov loops is used to generate a flux tube. In both cases 
the temporal extent has to be large enough to prevent strong contaminations by excited states 
respective finite temperature effects. What makes life difficult is the exponential fall-off of the 
signal with the temporal extent of loops. As already pointed out in the last section, in the dual 
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formulation in principle there is no problem to go to arbitrary low temperatures. We represent 
static charges by Polyakov loops and calculate the time averages of the interesting observables. 
The quality of the signal is constant for arbitrary time extents and for multiply charged systems, 
and also for increasing charge distances (as long as we use a lattice of same size). In this section 
we will investigate the flux tube profile and the total field energy for a system of (double) charges 
at a distance d. 

The connection between the flux tube cross section A and the string tension a was inves- 
tigated in ref. Assuming that flux tube models of confinement are consistent with lattice 
simulations, the following relation should be valid for £7(1) gauge theory with charges Q that 
are integer multiples of the coupling e: 

This means that if the cross section A is independent of the representation, the string tension 
will scale proportionally to the squared charge. Such a behaviour has been reported in ref. Q for 
three-dimensional lattice gauge theories. In a large class of phenomenological models, however, 
A is expected to increase with the representation. We are going to investigate the dependence 
of the flux tube cross section on the charge Q of the sources. For our analysis we will use the 
definition of A suggested in ref. |3| 

f&dA f E 2 dA 

where in the second expression we have exploited the exact validity of Gauss' law |5|, Let 
us first compare the flux tube profiles for single and double charges. In fig. [l| the longitudinal 
component E\\ of the electric field as well as its square E 2 are shown in the symmetry plane 
between charges as a function of transverse distance R from the axis. Note that the distance 
between the charges is 22 lattice spacings. In the upper plot it can be seen that the field strength 
of the double flux tube on the axis is significantly lower than twice the field strength of a single 
flux tube: The electric flux seems to spread over a wider transverse area. The plot for the 
distribution of E 2 looks very similar. We therefore expect that, calculating the cross sections 
according to definition (||), the double flux tube cross section will be significantly greater than 
that of a single flux tube, at least for this large charge distance. We show in fig. ^ that this is 
indeed the case. The ratio Adoubie/A is plotted as a function of the distance between charges. 
For larger distances it is equal to 2 within the errorbars; the deviation at small distances seems 
to result from Coulombic contributions important for these distances. 

In terms of the flux tube model (|?j) this increase of the cross section indicates that the 
string tension increases less rapidly than the squared charge. According to the ratio shown in 
fig. ^ we should expect a oc Q. We postpone a further discussion of this question to the next 
section, because the string tension should be extracted from the potential, which differs from 
the electromagnetic field energy. 

Nevertheless it is interesting to look at the total energy U of the electromagnetic field as 
a function of charge distance. This quantity is obtained by integrating the energy density 
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d = 22a 
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R [a] 

Figure 1: Flux tube profile of longitudinal electric field (E\\) and its square (£| ) in the symmetry 
plane between single (solid line) and double charges (dashed line) at a distance d = 22a for 
P = 0.96. 
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Figure 2: Ratio of cross sections of a double and a single flux tube in the symmetry plane as a 
function of charge distance d for (3 = 0.96. 

(E 2 + B 2 )/2 over the whole lattice^, i. e. summing up all plaquette contributions. The significant 
difference of the results between confinement and Coulomb phase is demonstrated in fig. [3|, where 
the total field energy U as a function of charge distance d is plotted for single and for double 
flux tubes. In the upper plot {(3 = 0.96) it can be seen that after a few lattice spacings the 
energy increases linearly, the ratio of slopes approximately equals 2. The lower plot shows the 
corresponding results at (3 = 1.10, where the system is in the Coulomb phase and may be 
approximated by classical electrodynamics. It is not possible to detect a linear increase, the 
ratio of energies is significantly larger (for (3 — > oo it is expected to approach 4). Simulations for 
this value of (3 have only been performed on smaller lattices, because quantum fluctuations are 
much larger and the energy nearly saturates at d ~ 5a. The ratio of the energies of a double and 
a single flux tube is shown in fig. again as a function of the charge distance d, for (3 = 0.96 
and (3 = 1.10. 

4 The energy of periodically closed flux tubes 

As already pointed out, by studying closed flux tubes one gets rid of end effects which makes it 
easier to determine the string tension. Instead of measuring the slope of the potential at large 
distances, we just have to measure the energy of the flux tube relative to the vacuum^. As we 

2 We use the Minkowski space notation here. The magnetic energy density relative to the vacuum is negative, 
therefore it comes to cancellations between electric and magnetic energy. 

^Strictly spoken the energy per unit length is not equal to the string tension because of string fluctuation 
contributions Jl3; this is not relevant, however, for the following discussion. 



8 



(3 = 0.96 




5 10 15 20 

d[a] 




Figure 3: Total field energy of a pair of charges (solid line) and a pair of double charges (dashed 
line) as a function of charge distance d for (3 = 0.96 and (3 = 1.10. 
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Figure 4: Ratio of the total field energies Udouble/U as a function of charge distance d for 
(3 = 0.96 and (3 = 1.10. In the Coulomb phase {(3 = 1.10) this ratio is difficult to determine for 
larger lattice sizes. 



will show this is a good quantity to distinguish between confined and deconfined phase. Further 
we want to work out the difference between the total field energy as calculated in the last section 
and the free energy calculated via a numerical integration over (3. 

In fig. H free energy and field energy per length are contrasted both for the three- and for the 
four-dimensional theory as a function of (3. A large difference between the quantities is observed 
at strong coupling. For (3 — ► the free energy diverges logarithmically (in both three and four 
dimensions), whereas the field energy per length tends to 1 (in units of the inverse squared 
lattice spacing). In the classical limit, however, both energies show the same behaviour. 

In three dimensions the free energy becomes smaller than the field energy at (3 ~ 0.9, as can 
be seen from fig. [|. In a range of weaker couplings where already many simulations have been 
performed in previous times, the discrepancy is approximately a factor 2 which is in agreement 
with ref. §|. In the four-dimensional theory the free energy is greater than the field energy for 
the whole confinement phase. At the phase transition both of them show a drastic decrease. 
On a lattice of finite transverse extent the energy is not exactly zero above the phase transition 
due to the classical energy of the corresponding homogeneous electric field on a torus. This 
kind of finite size effect may play a more important role in three dimensions as will be discussed 
in the next section. Thus we have seen that it is necessary to make a difference between free 
energy and total field energy, nevertheless these two lattice quantities can be related to each 
other by sum rules [11] which may be used to extract further informations on the system [12|. 
Corresponding calculations are in progress and will be reported elsewhere. 

Now we proceed to the investigation of closed double flux tubes in order to continue the 
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Figure 5: Free energy (solid line, obtained by numerical integration over f3) and total field energy 
(points) per unit length of periodically closed flux tubes as a function of the coupling (3. The 
left plot is for three-dimensional U(l) on a 8 2 x 16 lattice, the right plot for the four-dimensional 
theory on a 8 3 x 16 lattice. Errorbars are not drawn, since the estimated relative error for the 
free energy does not exceed a few percent, even in the critical /3-region of the phase transition. 



discussion of the last section. For the moment we focus on the four-dimensional theory; dif- 
ferences to the three-dimensional case will be discussed later. The free energy of a double flux 
tube shows a very similar /^-dependence as for the single flux tube. We are mostly interested 
in the ratio F^oubie/^] it is plotted in fig. |as a function of the coupling (3. In the confinement 
phase (starting from very small /3-values not shown in the plot) the ratio is 2 within the esti- 
mated errors, which could have been expected from the results of the last section. It can be also 
noticed, however, that between (3 ~ 0.9 and the phase transition at (3 ~ 1.01 this ratio drops 
significantly below 2. Above the phase transition, in the Coulomb phase, the ratio jumps to 4 
which is the classically expected result for the energy of a homogeneous field of double strength. 
The results shown in fig. |6| are for a 8 3 x 16 lattice; it has also been checked that there does not 
change much when considering larger spatial extents except that the statistics gets worse. 

How to understand these results in terms of the dual formulation? That the energy of a 
double flux tube is just twice the energy of a single flux tube in the strong coupling limit has 
a simple explanation. Let us start from a configuration with the electric flux identical to the 
double dual Dirac sheet (*n = 2 =^ *k = 2). The action of such doubly occupied plaquettes 
is higher than that of two singly occupied plaquettes for any finite value of (3. Therefore it is 
favourable for the system to form two completely separated flux strings (this is accomplished by 
a change in */) which in the ground state of the dual theory do not interact at allQ. For short 
open strings between doubly charged static sources, end effects cause a deviation from the ratio 
two. 

In ref. 0] it was argued that one should expect that the string tension scales like the squared 

4 It is not correct, however, to conclude from this "classical" behaviour of the dual theory that the confinement 
vacuum behaves like a dual type-II superconductor. As is shown in fig. ^| and will be analyzed later, the quantum 
fluctuations of the dual theory can lead to an attraction between the (7(1) flux tubes. 
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Figure 6: Ratio of free energies of doubly and singly charged closed flux tubes as a function of 
P for four-dimensional U(l) on a 8 3 x 16 lattice. 



charge through an extension of the Villain approximation used in ref. [p^[ , including multiply 
charged Wilson loops. This argument is only true, however, if the flux strings cover the same 
plaquette. As we explained, such configurations are unstable in the dual simulation. An analysis 
of the results obtained in ref. 0] for the three-dimensional theory will be performed in the next 
section. 

On the other hand, the classical Coulomb energy of a doubly charged state is four times 
higher. This can also be explained within the dual formulation. In the Coulomb phase the 
probability for higher plaquette values *k increases; for large (3 the /c-distribution may be ap- 
proximated by a continuous function p(k) which is symmetric with respect to k = if no charges 
are present. By introducing a closed single flux state on the lattice, this distribution is shifted 
a little by 8ko so that p(k) — > p(k — Sko); two flux quanta cause a shift by 2bk§. The opera- 
tor for the squared field strength introduced in eq. (|5|) in the weak coupling limit behaves as 
k 2 /2(3 + const (P). In this case the total energy of a flux state with respect to the vacuum is 
proportional to 

J dk [p(k - 5k ) - p(k)} g = fdk p{k) 

i. e. it scales with the squared charge of the system. 

A physically very interesting observation is the relative decrease of the double flux tube 
energy in the region of P approaching the phase transition. This indicates that there should be 
an attractive force between flux tubes in four dimensions. A closer investigation and discussion 
of this effect will be performed in section |(| 

Finally we also show the above discussed ratio between double and single flux tubes for the 
total field energy in fig. [?]. It qualitatively behaves very similarly to the ratio of free energies 
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Figure 7: Ratio of total field energies of doubly and singly charged closed flux tubes as a 
function of (3 for four-dimensional U(l). 



and also confirms the observed decrease below 2 in the confinement phase in the vicinity of the 
phase transition. Because the field energy of a flux state in the Coulomb phase is very small as 
discussed above, the errorbars for the corresponding ratios get very large. 

5 Analysis of the three-dimensional theory 

As already mentioned our results differ from those of ref. 0] which were obtained in a simulation 
of three-dimensional U(l) lattice gauge theory. Singly and doubly charged Wilson loops were 
calculated at = 2.4 on a 32 3 lattice, and the potential was found to scale with the squared 
charge within a few percent. Due to the nature of a standard calculation of correlation functions, 
the largest considered spatial extent of the loop was 6a. 

We will argue that the observed behaviour may be due to two effects which are relevant 
mainly for the three-dimensional theory: the finite transverse extent of the lattice (there is 
only one transverse direction for flux tubes in three dimensions) and the influence of the static 
sources, yielding a non-negligible Coulomb term in the potential (which is logarithmic in three 
dimensions) for small and intermediate charge distances. 

Let us first look again at the ratio of free energies for closed flux tubes. Since in three- 



dimensional U(l) the string tension is non-zero for all values of the coupling 14] we might expect 
that the ratio Fd ou bi e /F equals 2 for all values of (3. But if j3 is sufficiently large, the thickness 
of the flux tube approaches the transverse lattice size, and the cross section cannot increase 
proportionally to the charge any more. This causes a smooth transition of Fdoubie/ 'F from 2 to 
4. The dependence of this transition on the transverse extent N x of the lattice is illustrated in 
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Figure 8: The ratio of the free energies of doubly and singly charged closed flux tubes as a 
function of j3 in three-dimensional U(l) for various transverse lattice sizes N x . 

fig. ||. By enlarging N x the transition can be shifted to higher values of j3. Consider now (3 = 2.4 
and JVaj = 32 for comparison to ref. Q]: We are just in the transition region, the free energy of 
the closed double flux tube is approximately three times higher than that of the single flux tube. 
This shows that finite size effects already play a considerable role. We have also investigated 
the ratio of total field energies; again this ratio behaves very similarly. 

The Coulomb potential from static sources may further hide the proportionality of the string 
tension to the charge. Thus we have also computed the total field energy U of a system with 
static sources on a 32 3 lattice at = 2.4 and show the ratio Udoubie/U as a function of the 
charge distance in fig. ||. As demonstrated in the last section, we may expect that this ratio 
shows roughly the same behaviour as the one for the free energy. It can be seen that Udoubie/U 
indeed is compatible with 4 (but not with 2) at smaller distances. For larger distances it takes 
values in an intermediate region around 3, which is in agreement with the behaviour of closed 
flux tubes at the given parameters shown in fig. |8[ 

Thus we have shown that the data obtained by our dual simulations quite agree with data 
from standard simulations, but within the dual formulation it is much easier to exclude system- 
atic errors hiding the true nature of flux tubes. 

6 The interaction between flux tubes in four-dimensional U(l) 

For the three-dimensional theory our data have clearly demonstrated that the string tension 
scales proportionally to the charge, there is no detectable interaction between flux tubes. As 
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Figure 9: Ratio of total field energies for static double and single charge pairs as a function of 
distance d at j3 = 2.4 on a 32 3 lattice. 



already mentioned, the situation is different in four-dimensional U(l). In figs. |6| and ^ we have 
seen a small deviation from the proportionality to the charge in the /3-range slightly below the 
phase transition. We will now perform a closer investigation of the interaction between flux tubes 
for (3 = 0.96. We will try to clarify the origin of this interaction, to calculate the interaction 
energy in the case of flux tubes between static double charges, and to investigate its influence 
on the flux tube profiles. 

Let us first analyze the effect for periodically closed flux tubes. If we consider the total 



field energy, we are able to separate the interaction energy L 1 double — 217, 



single 



into individual 



components - electric and magnetic, longitudinal and transverse. The result is shown in fig. 10 
There is almost no interaction term for the longitudinal electric field -E| which gives the main 
contribution to the string tension. The transverse components of electric and magnetic field 
cancel each other, due to the symmetry of the lattice. As can be seen the main contribution of 
the interaction between flux tubes is due to the longitudinal mag netic field Bh This increases 
the average action but decreases the average energy of a double flux tube. As discussed in 
ref. H , the longitudinal magnetic field only plays an important role at higher values of j3 where 
the fluctuations in the dual model significantly increase. In terms of the dual simulation one can 
argue that these excitations of the string are enhanced by the presence of a second flux tube: 
It may happen that instead of 4 new plaquettes in a cube-like excitation only 2 plaquettes get 
occupied. It is plausible that such an opening of new channels also lowers the free energy of the 
system. The responsibility of B^ for the attractive interaction between flux tubes also makes 
clear why there is no such effect in three dimensions where we can only distinguish between 
three different types of fields (Eh, E± and B±). 
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Figure 10: The individual contributions to the interaction energy Udoubie — 2U S i ng i e f° r a peri- 
odically closed double flux tube at j3 = 0.96. The contribution of the electric field is positive, 
that of the magnetic field negative; therefore the total interaction energy is negative. 



If static double charges are attached to the ends of such a double flux tube, we expect an 
additional interaction term in the field energy, due to the repulsion between equal charges. The 



interaction energy Udoubie — 2U s i ng [ e is shown as a function of the flux tube length d in fig. |11 



A strong increase is observed till d = 2a, where Coulombic behaviour is dominating, for larger 
distances the interaction energy seems to decrease slightly (although difficult to resolve even 
in the dual simulation, this effect is in agreement with the observations for closed flux tubes). 
Thus the positive interaction energy for a doubly charged system is only due to the self energy 
of the double charges, for very long flux tubes we might expect that this is compensated by the 
negative contribution of the flux tube attraction. 

Finally we return to the investigation of flux tube profiles. Since two flux tubes were shown 
to attract each other, we should expect that two (initially) separated flux tubes prefer to form 
one single flux tube. To test this scenario we consider two charge pairs at a transverse distance 
of 4 lattice spacings in cc-direction. It is important to investigate long flux tubes which have 
sufficient mobility to "feel" the interaction. Thus we take again d = 22a and measure the 
distribution of the longitudinal electric field along the x-axis in the symmetry plane. The result 
is depicted in fig. 12. It differs significantly from the field which is obtained by a superposition of 
the contributions of two single flux tubes, shifted in transverse direction by ±2 lattice spacings. 
This demonstrates that the system seems to prefer the formation of one flux tube, in analogy 
to a dual type-I superconductor. 
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Figure 11: The interaction (field) energy of a pair of double charges as a function of charge 
distance at (3 = 0.96. 




Figure 12: Longitudinal electric field profile of two interacting flux tubes in the symmetry plane 
(E\\ , solid line). The length of flux tubes is d = 22a, the transverse distance of equal charges is 4a. 
For comparison, the dotted lines show the results for single flux tubes at x = —2a and x = +2a, 
and the dashed line corresponds to the superposition Ep_ + E\\ 2 of these two non-interacting flux 
tubes. 
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7 Conclusions 



The dual formulation of £7(1) lattice gauge theory has proved to be a very efficient tool for the 
numerical investigation of doubly charged systems and the interaction between flux tubes. We 
are able to exclude strong finite temperature effects by considering Polyakov loops with large time 
extents; long flux tubes become numerically accessible; doubly charged systems can be studied 
with equal accuracy; a simple implementation of periodically closed flux tubes eliminates the 
contributions from static sources; finally we can calculate the total field energy and the free 
energy (potential) of flux tubes by dual simulations. 

We have found that both in four and in three dimensions the string tension scales propor- 
tionally to the charge at strong coupling. Deviations from this behaviour in the energy of doubly 
charged systems in three-dimensional £7(1) have been shown to be due to finite size effects and 
to the influence of static sources. In four dimensions there is an attractive interaction between 
flux tubes, which becomes stronger towards the phase transition (a similar phenomenon has 
also been observed for Z2 gauge theory at criticality |0). This effect in £7(1) is due to the 
longitudinal component of the magnetic field; it can be "hidden" by the strong repulsion of the 
equal charges. For flux tubes of sufficient length we also demonstrated that a fusion into one 
flux tube is preferred. 

In ref. Q the result that the potential scales like the squared charge in three-dimensional 
£7(1) for the considered parameters was interpreted as a similarity between the confinement 
mechanism in Abelian and non-Abelian gauge theories. One should be very careful, however, 
in comparing the results. In SU(3) lattice gauge theory it was shown that the Casimir scaling 
of the string tension can only be observed up to intermediate distances, where the flux tube 
cannot change its character |15[ |. As the interquark distance increases, it becomes energetically 
favourable to pair-create gluons which screen higher representation sources. In addition, the 
phenomenon of flux tube coalescence, which is observed in non-Abelian gauge theories [16], is 
absent in U(l). 

Finally we want to comment on the interpretation of £7(1) lattice gauge theory in the con- 
fining phase as a dual superconductor. The numerical results we have presented show that flux 
tubes in four dimensions attract each other. This corresponds to the behaviour of a type-I 
superconductor. But this does not mean that £7 (1) is a dual superconductor in the sense that 
it is described exactly by a dual Higgs model. It is well-known that the £7(1) partition func- 
tion can be rewritten as a double limit (two couplings sent to infinity) of a dual Higgs model; 
as discussed in ref. [|| the corresponding transformation of expectation values shows that the 
agreement between the observables in the two models is not complete. Nevertheless, one can 
try to describe the observables of £7(1) lattice gauge theory approximately by an effective Higgs 
model with appropriately chosen parameters. According to our £7(1) results, we expect that 
such parameters should be those of a type-I superconductor. 
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